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Effects of Massive Neutrinos on the Large-Scale Structure of the 
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ABSTRACT 

Cosmological neutrinos strongly affect the evolution of the largest structures in the Universe, 
i.e. galaxies and galaxy clusters. We use large box-size full hydrodynamic simulations to in- 
vestigate the non-linear effects that massive neutrinos have on the spatial properties of cold 
dark matter (CDM) haloes. We quantify the difference with respect to the concordance ACDM 
model of the halo mass function and of the halo two-point correlation function. We model the 
redshift-space distortions and compute the errors on the linear distortion parameter p intro- 
duced if cosmological neutrinos are assumed to be massless. We find that, if not taken cor- 
rectly into account and depending on the total neutrino mass M v , these effects could lead to 
a potentially fake signature of modified gravity. Future nearly all-sky spectroscopic galaxy 
surveys will be able to constrain the neutrino mass if M v > 0.6 eV, using p measurements 
alone and independently of the value of the matter power spectrum normalisation Og. In com- 
bination with other cosmological probes, this will strengthen neutrino mass constraints and 
help breaking parameter degeneracies. 
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1 INTRODUCTION 

Neutrinos are so far the only dark matter candidates that we actu- 
ally know to exist. Since deviations from the stand ard model in the 
form of extra neutrino species are st ill uncertain dGiusarma et al.l 
1201 ll ; iGonzalez-Morales et alj|201lh . in this paper we focus on 
standard neutrino families only. It is now established from solar, 
atmospheric, reactor and accelerator neutrino experiments that neu- 
trinos have non-zero mass implying a lower limit on the total neu- 
trino mass given by M v = £ra v ~ 0.05 eV (Lesgourgues & Pastor 
2006), where m v is the mass of a single neutrino species. On 
the other hand the absolute masses are still unknown. Since neu- 
trino mass affects the evolution of the Universe in several ob- 
servable ways, its measurements can be obtained from different 
cosmological probes as observations of Cosmic Microwave Back- 
groun d (CMB), galaxy clu stering, Lyoc forest, and weak lensing 
data dAbazaiian et al.ll20lTh . 

In particular, a thermal neutrino relic component in the Uni- 
verse impacts both the expansion history and the growth of cosmic 
structures. Neutrinos with mass < 0.6 eV become non-relativistic 
after the epoch of recombination probed by the CMB, and this 



mechanism allows massive n eutrinos to alter the matter-ra diation 
equality for a fixed Q.,„h 2 dLesgourgues & Pastor] 120061) . Mas- 
sive neutrinos act as non-relativistic particles on scales k > k m = 
0.018(m v /leV) 1 / 2 ii, 1 ,/ 2 /zMpc" 1 , where k m is the wave-number 
corresponding to the Hubble horizon size at the epoch z m when the 
given neutrino species becomes non-relativistic, £l m is the matter 
energy density and h = Z/o/lOOkms^'Mpc -1 . The large velocity 
dispersion of non-relativistic neutrinos suppresses the formation of 
neutrino perturbations in a way that depends on m v and redshift z, 
leaving an imprint on the matter power spectrum fo r scales k > 
O.S2H(z)/H n /(l + z) 2 (m v /leV) /iMpc' 1 jTakada et al.l 



l2006tlLesgourgues & Pastorl2006l) . where neutrinos cannot cluster 
and do not contribute to the gravitational potential wells produced 
by cold dark matter and baryons. This modifies the shape of the 



(see e.g.|5proshkevich et al. 1981 : Hu et al.lll998l: 


Abazajian et al. 


2005; Kiakotou et al. 2008; Brandbyge et al. 2010; 


Viel et al. 2010, 



and reference therein). 

Massive neutrinos affect also the CMB statisti cs. WMAP7 
alone constrains My < 1.3 eV (Komatsu et al. 2009) and, thanks 
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to the improved sensitivity to polarisation and to the angular power 
spectrum damping tail, forecasts for the Planck satellite alone give 
a l-o" error on the total neutrino mass of ~ 0.2 — 0.4 eV, de- 
pending o n the assumed cosmological model and fiducial neutrino 
mass (e.g. iPerotto et al.ll20ol : iKitching et alj|2008t and references 
therein). Moreover, the combination of present data-sets from CMB 
and large-scale structure (LSS) yields an upper limit of M v < 0.3 
eV dWang et alj|2005k l^khlinin et alj|2009l; iThomas et alj|20ld : 
iGonzalez-Garcia et alj|201(t iReid et alJkoid. e.g.). A further ro- 
bust constraint on neutrino masses has been obtained using the 
Sloan Digital Sky Survey flux powe r spectrum alone , finding an up- 
per limit of M v < 0.9 eV (2a C. L.) dViel et al.l2oTol) . However, the 
tightest contraints to date in terms of a 2a upper limit on the neu- 
trino masses have been obtained by combining the Sloan Digital 
Sky Survey flux power from the Lyman alpha forest wi th CMB and 
alaxy clustering data and result in £m v < 0.17 eV 1 Selia ket al.l 
20061) . Somewhat less constraining bounds have been obtained by 



Goobar et al.l ( l2006T) . while for forec asting on future joint CMB and 



Lyman-alpha constraints we refer to lGratton et al.ld2008l) . For fur- 
ther dis cussion on neutr i no ma ss constraints from different probes 
see e.g. lAbazaiian et alj ( 1201 II) and reference therein. 

The forecasted sensitivity of future LSS experiments, when 
combined with Planck CMB priors, indicates that observations 
should soon be able to detect signatures of the cosmic neutrino 
background and measure the neutrino mass even in the case of 
the minimum mass M v 
Kitching et all 120081; 



: 0.05 eV (e.g. Hannestad & Won 



LSST Science Collaborations et al 



2007 



2009 



Hannestad l20ld : lLahav et alj|20ict) . In particular. ICarbone et all 
| |2011 ) show that future spectroscopic galaxy surveys, such as EU- 
CLID, JEDI and WFIRST, not only will be able to measure the 
dark-energy equation of state with high accuracy, but they will de- 
termine the neutrino mass scale independently of flatness assump- 
tions and dark energy parametrization, if the total neutrino mass 
My is > 0. 1 eV. On the other hand, if M v is < 0. 1 eV, the sum of 
neutrino masses, and in particular the minimum neutrino mass re- 
quired by neutrino oscillations, can be measured in the context of a 
ACDM model. 

It is therefore mandatory to measure with high accuracy the 
growth history of large-scale structures in order to obtain the nec- 
essary cosmological information, excluding possible systematics 
due to the incorrect assumption that neutrinos are massless. One 
way of determining the growth of structure is through the red shift- 
space distortions (RSD) of the galaxy distribution, caused by the 
line-of-sight component of galaxy peculiar velocities. RSD can be 
exploited in large deep redshift surveys to measure (if the galaxy 
bias is measured independently) the growth rate of density fluc- 
tuations / = dlnD/dlna, with D being the linear density growth 
factor and a = 1/(1 + z), or to mea sure the linear redshift-space 
distortion parameter (3 dKaisedfl987l) that depends on the growth 
rate / and the galaxy linear bias b. In particular, (3 in the presence 
of massive neutrinos depends on both redshift and wave-numbers 
(3(z,fe) = f(z,k)/b(z), since in this case the linear growth rate 
f(z,k), being suppressed by free-streaming neutrinos, acquires a 
scale dependence already at the linear level (Kia kotou et al. 1120081) . 

There are two types of RSD with competing effects acting 
along opposite directions on the observed galaxy correlation func- 
tion. While, for large separations, large-scale bulk peculiar veloc- 
ities produce a flattening effect on the correlation function and 
give information on the growth of structures, on small scales, ran- 
dom peculiar velocities cause the so-called Fingers o f God (FoG), 
stretching compact structure s along the line-of-sight dScoccimarrd 
l2004ISong & Percivalll2009l) . 



RSD have bee n the subject of many analyses, as reviewed 
in Hamilton ll998l) . The latest large galaxy surveys that have 
enabled measurements of RSD via the correlation function and 
the po wer spectrum are th e 2-degree Field Galaxy Redshift 



Survey jPeacock et alj l200ll ; lHawkins et al 



20041) and the Slo a n Digital Sky Survey 1 Tegmark et al.l 
Zehavi et alj 120051; iTegmark et all [200 3; lokumura et al.l 



20031; IPercival et al.l 



2004 



2008 



Cabre & Gaztanagal l2009al |b|). More over, also the VIM OS-VLT 



Deep Survey have been exploited in lGuzzo et al. (2008) for RSD 
determinations from the correlation function. 

Since the linear theory description is valid only at very large 
scales, an extension of the theoretical description has been at- 
tempted to quasi-linear and non-linear scales u sing empirical meth- 
ods based on the so-called streaming model (Peebles 1980), con- 
sisting of linear theory and a convolution on the line-of-sight with 
a velocity distribution. This model describes the FoG elongation 
along the line-of-sig ht due to random moti ons of viriali s ed ob - 
jects jjacksonlll972l). It ha s bee n shown by | Guzzo et ail d2008l) . 
ICabre & Gaztafiagal J2009al) and IPercival & Whitd d2009b that on 
quasi-linear scales a streaming model with a Gaussian velocity dis- 
persion is a good general fit to the redshift-space power spectrum. 
Howev er, this model is not accurate on very small and very large 
scales dTaruva et al.ll2010T ; IOkumura & Jinj|201ll : lRaccanelli et al.l 
l20ld) and fitting functions based on simulation results have been 
used faatton & Coldll999l;IScoccimarrdl2004lTinker et al.ll200d: 
Tinkerl l2007l : iToc chini-Val entini et alj|201ll) . Anyway, to the pur- 
pose of this paper, the streaming model is accurate enough to ro- 
bustly constrain the effect of massive neutrinos on RSD when ap- 
plied on scales < 50 h~ 1 Mpc . 

In this work, we compare analytic results against a set of large 
N-body hydrodynamical simulations developed with an extended 
version o f GADGET III, th at is an improved version of the code de- 
scribed in Springelld2005t) . further modified to take into account the 
effect of massi ve free-streaming neutrinos on the evolution of cos- 
mic structures dViel etalJ boiO). It is well known that galaxy/halo 
bias on large, linear scales is scale-independent, but becomes non- 
linear and therefore scale-dependent on smaller scales. This effect 
can be mimicked or enhanced by the presence of massive neutrinos. 
Therefore, the effect of massive neutrinos on the galaxy clustering 
in the quasi non-linear regime has to be explored via N-body sim- 
ulations to encompass all the relevant effects, and analyse possible 
sources of systematic errors due to non-linearities and galaxy bias 
scale-dependence. In particular, in this work we will focus on the 
DM halo mass function (MF), the DM halo bias and RSD. 

The rest of the paper is organised as follows. In fj2]we review 
our method and the adopted modelling of RSD. In ij3]we describe 
the exploited set of N-body simulations and present our results on 
the neutrino effects on LSS. Finally in *]5]we draw our conclusions. 



2 FORMALISM TO MODEL REDSHIFT DISTORTIONS 
2.1 Overview 

In this section we describe how RDS are generated in the observed 
galaxy correlation function. An observed galaxy redshift is com- 
posed by the two additive terms, 



£obs — £c " 



VII 

c 



(1) 



where z c is the cosmological redshift, due to the Hubble flow. The 
second term of Eq. |QJ is caused by galaxy peculiar velocities where 
v || is the component parallel to the line-of-sight. The real comoving 
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distance of a galaxy is given by 



"II 



Jo H( z ' c y 



(2) 



where H(z' c ) is the so-called Hubble rate. When the distances are 
computed replacing z c with z a b s m Eq. (O, i.e. without correcting 
for the peculiar velocity contribution, we say to be in the redshift- 
space. We will refer to the redshift-space spatial coordinates using 
the vector s, while we will use r to indicate the real-space coordi- 
nates. 

Fundamental information is hinted in the anisotropies of an 
observed galaxy map in redshift-space. A useful statistics widely 
used to describe the spatial properties of a general astronomical 
population is the two-point correlation function, \{r), implicitly de- 
fined as dP\2 = n 2 [l +!^(r)]dV\dV2, where dP\2 is the probability 
of finding a pair with one object in the volume dV\ and the other 
in the volume dVi, separated by a comoving distance r. It is con- 
venient to decompose the distances into the two components per- 
pendicular and parallel to the line-of-sight, r = (r^,r^), so that the 
correlation becomes a two-dimensional function of these variables. 
When measured in real-space, the contour lines of |(rj_ , rii ) are cir- 
cles for an isotropic population of objects as galaxies. Instead, in 
redshift-space % is distorted: at small scales (<1 h Mpc ) the dis- 
tortion is caused by the random motions of galaxies moving inside 
virialised structures. This motion changes the shape of t, in the di- 
rection parallel to the line-of-sight, producing the observed FoG. 
At large scales the coherent bulk motion of virialising structures 
squashes the correlation function t, perpendicularly to the line-of- 
sight. 

A different kind of dist ortion, called geometrical o r Alcock- 
Paczynski (AP) distortion dAlcock & Paczvnskil 1 19791) , can be 
present if the cosmological parameters assumed in Eq. (|2} are not 
the same as the true cosmological model of the Universe. In what 
follows, we will not consider this effect since massive neutrinos do 
not produce a geometrical distortion if the present-day total matter 
energy density parameter £2 m is held fixed. 



2.2 Modelling the dynamical distortions 

At large scales and in the plane-parallel approximation, the dynam- 
ical distortions can be parameterised in the Fourier space as follows 



p(k) = (i + fr 2 ) 2 p lin (k) 



(3) 



where P\\ n (k) is the linear power spectrum of the matter density 
fluctuations and jj is the cosine of the angle between k and the line- 
of-sight. Fourier transforming equation ((3] gives 

%{s,n) = ^(s)P (v) + ^)P2M + ^(s)P4(v), (4) 

where the func tions P; represent the Legendre polynomials 
(Hamilton 1992). The multipoles t, n (s), n = 0,2,4, can be written 
as follows 



$4W 



8p^ 
35 



%{r) + \%{r)-\%{r) 



(5) 



(6) 



(7) 



where t,(r) is the real-space correlation function, and the barred 
correlation functions are defined as 



%{r) = 4 fdr'^r'y 1 
r 3 Jo 

f(r) = ^ [dr'Wy 4 



(8) 



(9) 



Eq. 10 can be used to approximate the correlation function at large 
scales. To include in the model also the small scales, as discussed 
in SjT] we adopt the streaming model and use the following formula 



S(*x,*||) = f^f{ V )%{s ± ,s {l - V /H(z)/a(z)), 



(10) 



where /(v) is the distribution function of random pairwise veloci- 
ties that are measured in physical (not comovi ng) coordinates (but 
see e.g. IScoccimarrol 12004: Matsubaral l2004h . In this work, we 
adopt for /(v) the form 

V2\v\ 



/(v) 



1 



: exp 



rJi2\/2 * \ CJ|; 

where o"i2 is the dispersion in the pairwise peculiar velocities. 



(11) 



3 THE ADOPTED SET OF NUMERICAL SIMULATIONS 

The set of simulation s we c onsider in this work have been per- 
formed by IViel et al.l ( 1201 Of) with the hydrodynamical TreePM- 
SPH (Tree Particle Mesh-Smoothed Particle Hydrodynamics) code 
GADGET in, whic h is an improv ed and extended version of the 
code described in ISpringell d2005l) . This code has been modified 
in order to simulate the evolution of the neutrino density distribu- 
tion. The cosmological model adopted in the simulations is based 
on cold dark matter and assumes the presence of the cosmolog- 
ical constant (ACDM): n s = 1, Q m = 0.3, Q. b = 0.05, £2 A = 0.7 
and h = 0.7 (Hq = 100/z km/s), plus a cosmological massive neu- 
trino component fl v = M v /(/z 2 93.8eV) (ACDM+v). In what fol- 
lows, we consider only the so-cal led "grid b ased implementation" 
of the simulations developed by IViel et ai] d2010h. where neutri- 
nos a re treated as a fluid (see also Brandbyge & Hannestadl l2009l 
l2010f) . In this implementation the linear growth of the perturbations 
in the neutrino component is followed by interfacing the hydro- 
dynamical code w ith the public available Boltzmann code CAME0 
dLewis et alj2000l) . More specifically, the power spectra of the neu- 
trino density component are interpolated in a table produced via 
CAMB of one hundred redshifts in total, spanning logarithmically 
the range z = — 49. The gravitational potential is calculated at the 
mesh points and the neutrino contribution is added when forces are 
calculated by differentiating this potential. 

In this approach the gravitational force due to neutrinos is cal- 
culated based on the linearly evolved density distribution of the 
neutrinos in Fourier space. This implementation has the advantage 
that it does not suffer from significant shot noise on small scales, 
yielding therefore higher accuracy at scales and redshifts where the 
effect of the non-linear neutrino evolution is still moderate, espe- 
cially for small neutrino masses. Further advantages of such a grid 
based approach, aside from eliminating the Poisson noise, are the 
reduced requirements with regard to memory (there are no neutrino 
positions and velocities to be stored) and computational time. 

The initial conditions of this set of simulations were generated 
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Figure 1. DM halo MF as a function of M v and redshift. Left: MF of the SUBFIND haloes in the ACDM N-body simulation (blue circles) and in the 
two simulations with M v = 0.3 eV (magenta triangles) and M v = 0.6 eV (red squares). The blue, magenta and red lines show the halo MF predicted by 
ISheth & Tormenl 1*2002), where the variance in the density fluctuation field, o"(M), for the three cases, M v = 0,0.3,0.6 eV, has been computed using the linear 
matter P(k) extracted from CAMB. Right: ratio between the halo MFs of the simulations with and without neutrinos. The green triangles show the MF ratios 
of the FoF haloes, while yellow circles show the ones of the SUBFIND haloes. The lines represent the ST-MF ratios: the black solid lines are the MF ratios 
predicted for M v =0.3 eV and M v = 0.6 eV; the red dashed lines are the same ratios but assuming p = p c ■ (Q m - fly) in the ST-MF formula Eq. CTD 
IBrandbvge et al J20ld) : finally, the blue dotted lines are the ratios between the ST MFs in two ACDM cosmologies, which differ for the 0g normalisation, as 
explained in the text. The error bars represent the statistical Poisson noise. 



based on linear matter power spectra separately computed for each 
component (dark matter, gas and neutrinos) with CAMB. The total 
matter power spectrum was normalised such that its amplitude (ex- 
pressed in terms of ffg) matches the CAMB prediction at the same 
redshift. The mass per simulation particle at our default resolution 
is 1 .4 • W l0 M e /h and 6.9 ■ 1O 1o M //i for gas and dark matter, re- 
spectively. In this work we consider the set of simulations with a 
box of comoving volume V = (512/i -1 Mpc) 3 , and total neutrino 
mass M v = 0,0.3,0.6 eV, respectively. 



To identify DM haloes and their substructures we have 
used two different algorithms: a standard Friends-of-Friends (FoF) 
group-finder with lin king length b = 0.2, and the SUBFIND al- 
gorithm described in ISpringel et all fcOOlh . Apart from the right 
panel of Fig. [T] all the results presented in this paper have been 
obtained using our sub-halo catalogues, composed by the gravita- 
tionally bound substructures that SUBFIND identifies in each FoF 
halos. However, as we have explicitly verified, all the main conclu- 
sions of this work do not change if we consid er the halo catalogue s 
instead. This happens because, as showed bv lGiocoli etail feOlCh . 
in the mass range considered in this work the total mass function 
of haloes and sub-haloes is manly dominated by the halo systems. 
At z — the sub-halo contribution start to be seen only for masses 
< 10 10 M© /h. For all the considered M v values, we have restricted 
our analysis in the mass range M m ; n < M < M max , where M m j n = 
2 • 1O 12 M //i and M max = 2 • 10 15 ,5 • 10 14 ,3 • 10 14 , 1O 14 M //i at 
z — 0,0.66, 1,2, respectively. 



4 RESULTS 

In this Section we show how the halo MF, the halo bias, and the lin- 
ear redshift-space distortion parameter P get modified with respect 
to the standard ACDM case, when a massive neutrino component 
is taken correctly into account. In particular we compare the results 
between the ACDM and the ACDM+v cosmologies, and analyse if 
our findings agree with analytical predictions in the literature. 



4.1 The halo mass function 

As mentioned in iJT] the free- streaming of non-relativistic neutri- 
nos contrasts the gravitational collapse which is the basis of cos- 
mic structure formation. The first consequence of this mechanism 
is represented by a significant suppression in the average number 
density of massive structures. This effect can be observed in the 
high mass tail of the halo MF as measured from our set of simula- 
tions, and shown by the data points in the left panel of Fig.[T] For 
a fixed amplitude of the primordial curvature perturbations A^, the 
amount of the number density suppression depends on the value of 
the total neutrino mass M v . From the comparison of the correspond- 
ing MFs, we recover what theoretically expected, i.e. the higher 
the neutrino mass is, the larger the suppression in the comoving 
number density of DM haloes becomes. The suppression affects 
mainly haloes of mass 1O 14 M //i < M < 10 15 M Q //i, depending 
slightl y on the redshift z - This result is in agreement with the find- 
ings of IBrandbvge et alj j2010l) . In the same plot , we compare the 
measu red MFs with the analytical predictions of She th & Tormenl 
d2002l) (ST), represented by the solid, dotted and dashed curves, 
corresponding to the values M v = 0,0.3,0.6 eV, respectively. The 
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Figure 2. Real-space two-point auto-correlation function of the DM haloes in the ACDM N-body simulation (blue circles) and in the simulation with M v = 0.6 
eV (red squares). The blue and red lines show the D M correlation funct ion, for M v = and M v = 0.6 eV, respe ctively, obtained b y Fourier trasforming the 
non-linear power spectrum extracted from CAMB (Lewis & Bridle 2002) which exploits the HALOFIT routine I Smith et al. 2003). The bottom panels show 
the ratio between the halo cor relation functi on of the simulations with and without neutrinos. The error bars represent the statistical Poisson noise corrected at 
large scales as prescribed bv lMo et al.Nl992h . 



ST fit is based on the f act that the halo MF can be written as 
jPress & SchechteJll974l) 



MdMdn{M,z) _ rf(n di ^ 
p dM 



(12) 



with £ = [8 sc (z)/o(M)] 2 , where § sc (z) = 1.686 is the overdensity 
required for spherical collapse at z in a ACDM cosmology, and p = 
£2;«Pc, where p c is the critical density of the Universe. Here £2 m = 
&cdm + Qb + &v, and dn(M,z) is the number density of haloes in 



the mass interval M to M + dM. The variance of the linear density 
field, <J 2 (M), is given by 



2 (M) = I dk 



k 2 P\in(k) 
2tz 2 



W(kR)\\ 



(13) 



where the top-hat window function is W(x) = (3/x 3 )(sinx — 
xcosx), with A = (3M/47IP) 1 / 3 . 
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Figure 3. DM halo bias, b = (cjhaio/^DM) ' 5 , measured from the ACDM 
simulations (blue circles) and from the two simulations with M v = 0.3 eV 
(magenta triangles) and M v = 0.6 eV (red squares). The error bars repre- 
sent the propaga ted Poisson noise corrected at large scales as p rescribed by 
iMo et al] fl992h . Dotted lines are the theoretical predictions of S heth et al] 
1200 ll) (Eg. {B}). 

The four panels show the results at different redshifts, as 

labeled. 



Figure 4. Mean bias (averaged in 10/i~' Mpc < r < 50 h~ [ Mpc) as a func- 
tion of redshift compared to the theoretical predictions of Sheth et al. 1 2001) 
(dotted lines) (Eq. (T3J). Here the dashed lines represent the theoretical ex- 
pectations for a ACDM cosmology renormalized with the Gg value of the 
simulations with a massive neutrino component. The error bars rep resent the 
propagated Poisson noise corrected at large scales as prescribed bv lMo et"al] 
U992). 



The ST fit to C/(0 is 

/ 1 \ /C'\ I/2 e -^'/2 
C/(C)=a(i + ^) (I) - 7r , (14) 

with = 0.707^, p = 0.3, and A = 0.3222 determined from the 
integral constraint / /(Qrf^ = 1. 

In Fig. Q] the variance of the density field, a 2 (M), has been 
computed with th e matter power spectrum extracted from CAMB 
dLewis et alj2000h . using the same cosmological parameters of the 
simulations. In particular, in the left panel we show the MF of 
sub-structures identified using the SUBFIND algorithm, where the 
normalisation of the matter power spectrum is fixed by A^(fcrj) = 

2.3 x 10~ 9 at k = 0.002/Mpc dLarson et al]|201 ll) . chosen to have 
the same value both in the ACDM+v and in the ACDM cosmolo- 
gies. The error bars represent the statistical Poisson noise. 

In the right panel of Fig. [T] we show the ratios between the 
halo MFs evaluated in the ACDM+v and ACDM cosmologies. In 
particular, the triangles represent the ratios of the halo MFs mea- 
sured directly from the simulations afte r identifying the structures 
via a FoF group finder (Springel 2005), and the filled circles rep- 
resent the ratios between the ACDM+v and ACDM MFs for the 
substructure evaluated via the SUBFIND algorithm. 

The curves in the panel show the corresponding ST predic- 
tions for three cases: the solid lines represent the ratios be- 
tween the theoretical MFs when the total Cl m = Q. c d m + £lb + 
is inserted in Eq. d 1 2b through the expression p = £2 m p c ; ii) the 
dashed lines represent the theoretical MF ratios when the quantity 
p = (0. c dm + &b)Pc is used in Eq. J 1 2b ; Hi) finally, the dotted lines 
represent the ratios between the ST MFs in two ACDM cosmolo- 
gies, which differ for the as normalisation, i.e. the ratio between 
the MF in a ACDM cosmology having the same as value of the 



simulations with a massive neutrino component, and the MF in a 
ACDM cosmology having a as value in agreement with the CMB 
normalisation A^(fco) = 2.3 x 10~ 9 and in the absence of massive 
neutrinos. 

We note that the MFs of the haloes obtained with the FoF al- 
gorithm look to be better fitted by the theoretical predictions of the 
;)-case, while the MFs of the substructures obtained with the SUB- 
FIND algorith m have a trend much mo re similar to the predictions 
of the !f)-case ( iBrandbyge etaitold) . Moreover, we see that, as 
the redshift increases, the suppression of the halo number density 
due massive neutrinos moves also towards masses M ^ 10 14 Mq / 'h. 
As an example, the number density of haloes with mass W u M Q /h 
at z = decreases by ~ 15% for M v = 0.3 eV and by ~ 30% for 
M v = 0.6 eV, and, at z = 1, b y~ 40% and ~ 70 %, respectively. 

As already discussed in IViel et all j20ich . the free-straming 
of massive neutrinos leads to the well known degeneracy between 
the values of as and M v . However, as the »()-case lines show, 
the difference between the MFs with and without neutrinos does 
not reduce merely to a as renormalisation of the background cos- 
mology, since, even renormalising to the same as, neutrinos free- 
streaming alters the MF, changing its shape and amplitude espe- 
cially for the less massive objects (compare the dotted-blue and 
solid-black lines in the right panel of Fig. [T). However, the possi- 
bility to measure this effect and break the M v -Cs degeneracy de- 
pends on the value of the neutrino mass and on the sensitivity in 
measuring the MF at masses M < 10 14 Mq//!. The resolution lim- 
its of our simulations do not allow us to study the properties of 
haloes with masses M < W l2 Ma,/h. Future galaxy surveys, like 
EUCLID, should be able to break this degeneracy measuring the 
number counts of low mass galaxies to gether with the overa ll shape 
of the matter power spectrum (see e.g. lcarbone et al.l201 lb . More- 
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over, even with present observational data it has al ready been possi- 
ble to break the M v -a% degeneracy. For instance, IViel et alJ l Hoiol) 
robustly constrained the neutrino masses using the Sloan Digital 
Sky Survey flux power spectrum alone, without CMB priors on Og. 



4.2 The halo clustering and bias 

As well known, massive neutrinos strongly affect also the spatial 
clustering of cosmic structures. As explained in 612.21 a standard 
statistics generally used to quantify the clustering degree of a pop- 
ulation of sources is the two-point auto-correlation function. Al- 
though the free-streaming of massive neutrinos causes a suppres- 
sion of the matter power spectrum on scales k larger than the neu- 
trino free-streaming scale kf s , the halo bias results to be signifi- 
cantly enhanced. This effect can be physically explained thinking 
that, starting from the same A^(feo) as initial condition, due to the 
suppression of massive neutrino perturbations, the same halo bias 
would correspond, in a ACDM cosmology without neutrinos, to 
more massive haloes, which, as well known, are typically more 
clustered. 

In fact, Fig. [2] shows, at different redshift s, the two-point 
DM h alo correlation function measured using the lLandv & S zalav 
dl993h estimator, compared to the correlation function of the matter 
density perturbations. We observe that, while for a fixed A^(&o), 
due to neutrino free-streaming, the total matter correlation func- 
tion decreases with respect to the ACDM case, especially on small 
scales (compare the solid-blue and dashed-red lines in Fig. [2}, the 
halo correlation function undergoes the opposite trend (compare 
the data points in Fig. O, so that the matter perturbation suppres- 
sion is in some way compensated by a stronger spatial clustering of 
the massive haloes. 

In particular, the halo clustering difference between the 
ACDM and ACDM+v cosmologies increases with the redshift (as 
it happens also for the halo MFs). For M v = 0.6 eV we find that 
the halo correlation function in the presence of massive neutrinos 
at z = 1 is ~ 20% larger than in a pure ACDM model, and at z — 2 
the difference rises up to ~ 40% (see the bottom panels of Fig.[2}. 

This effect is even more evident in Fig. [3] and [4] that show 
the effective bias measured from the simulations (symbols) com- 
par ed to the analytica l predictions (dotted lines), obtained using 
the lShethetall feOOlt) (SMT) bias, weighted with the ST MF of 
Eq. (12): 



Hz) 



f^n(M,z)b SMT (M,z)dM 



(15) 



Q Z n{M,z)dM 

where M m j n and M max have been defined in $3] Also in this case, the 
theoretical expectations reproduce correctly the numerical findings, 
inside the statistical errors, and, as in the ACDM cosmology, the 
halo bias results to be scale-independent on large scales, while the 
effect of non-linearities starts to be important for separations r < 20 
h^ 1 Mpc . 



4.3 Redshift-space distortions 

As it happens for the halo MFs and clustering, also RSD are 
strongly affected by free-streaming neutrinos. Fig.[5]shows the real 
and redshift-space correlation functions of DM haloes extracted 
from the simulations as a function of the neutrino mass. In the 
presence of massive neutrinos the rms of galaxy peculiar veloci- 
ties is smaller than in a pure ACDM cosmology, due to the sup- 
pression of both the growth rate f(k, z) and the matter power spec- 



trum P(k,z), which enter the bulk flow pred icted by linear theory 
dKiakotou et alJl2OO8l : lEigar0v & Lahavll2005l) : 

(v 2 (R*)) = (27T 2 )- 1 H 2 Jdkf 2 P hn (k)Wc(kR*) , (16) 

where Wc(kR*) is the window function, e.g., for a Gaussian sphere 
of radius R*, W(kR t ) = exp(—k 2 R 2 /2). This effect competes with 
the increase of the halo bias discussed in £14.21 resulting in a 
redshift-space halo correlation function slightly suppressed in a 
ACDM+v cosmology. In the bottom panels of Fig.[5]we show the 
ratios t,(s)/t,(r) compared to the theoretical values represented by 
the large-scale limit of Eq. (4) 

M-l + ? + f <>7) 

W) 3 5 

The effect of massive neutrinos on RSD is evident in particu- 
lar when the correlation function is measured as a function of the 
two directions perpendicular and parallel to the line-of-sight. In 
fact, from the top and bottom panels of Fig. [6] we observe that, 
in the case of massive neutrinos, the spatial halo clustering is less 
enhanced in redshift-space than in real-space. On large scales, this 
effect is due to the lower value of (v 2 (R*)} when neutrinos free- 
streaming is taken into account. On small scales, our analysis shows 
that also FoG get decreased in the presence of massive neutrinos, so 
that the best-fit values of (3 and O12, derived by modelling galaxy 
clustering anisotropies, result to be different than what expected 
in a ACDM cosmology. This might induce a bias in the inferred 
growth rate from data analysis, and therefore a potential ly false 
signature of modified gravity (see e.g. Simpson e t al.l201 H) . More- 
over, estimates of p and 012, when compared with the ACDM ex- 
pectations, yield an indirect neutrino mass measurement and may 
help breaking degeneracies with the other cosmological parame- 
ters. 

We quantify these effects in Fig. |7j which shows the best-fit 
values of p and O12 as a function of M v and z, where we have ne- 
glected their scale-dependence which, for the neutrino masses con- 
sidered in this work, is small enough that stat istical errors hide devi - 
ations of p and (J12 from spatial uniformity (Kia kotou et alj|2008h . 
Therefore, in this case we have considered the linear redshift- 
space distortion parameter as a function of the redshift alone, 
p = f(z)/b(z) ^ a m (z) y /b(z), with y= 0.545. The data points of 
Fig. UJ show that neutrinos free-streaming suppresses p and O12 
by an amount which increases with M v and z, and, fixed A^(&o), 
is clearly distinguishable from the corresponding ACDM values 
(dashed lines). As an example, at z = 0.6 the p best-fit values de- 
crease by ~ 10% for M v = 0.3 eV, and by ~ 25% for M v = 0.6 eV. 
Likewise, the O12 best-fit values decrease by ~ 25% for M v = 0.3 
eV, and by ~ 45% for M v = 0.6. 

On the other hand, the P best-fit values fall in the shaded grey 
bands, which represent the propagated ~ 10% theoretical bias er- 
ror. These bands contain also the theoretical predictions obtained in 
a ACDM cosmology renormalised with the rjg value of the simula- 
tions with a massive neutrino component (blue dotted lines). This 
means that, if an error of ~ 10% is assumed on bias measurements, 
we are not able to distinguish the effect of massive neutrinos on 
P when the two cosmological models with and without v are nor- 
malised to the same Og. 

In Fig.[8]we show, as a function of M v and z, the relative differ- 
ence between the theoretical p values calculated in the ACDM+v 
and ACDM cosmologies, normalised to the same Gg. At z = 1 and 
for M v > 0.6 eV, the relative difference with respect to the M v = 
case is AP/P > 3%. This result is interesting, since future spectro- 
scopic galaxy surveys, as EUCLID, JEDI and WFIRST, should be 
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Figure 5. Two-point auto-correlation function in real and redshift space of the DM haloes in the ACDM N-body simulation (blue circles) and in the simulation 
with My = 0.6 eV (red squares). The bottom panels show the ratio between them, compared with the theoretical expectation given by Eq. (TTJ. The error bars 
represent the statistical Poisson noise corrected at large scales as prescribed by Mo e t al.lfl992l) . 



able to measure the linear redshift-space distortion parameter with 
errors ^ 3% at z ^ 1, per redshift bin. 



5 CONCLUSIONS 

In this work we have studied the effect of cosmological neutrinos 
on the DM halo mass function, clustering properties and redshift- 
space distortions. To this purpose we have exploited the grid imple- 
me ntation of t he hy drodynamical N-body simulations developed 
by IViel et alj | |2010|) . which include a massive neutrino compo- 



nent, taking into account the effect of neutrinos free-streaming on 
the cosmic structure evolution. In ord er to model RS D, we have 
adopted the so-called streaming model (Peebles 1980), which con- 
sists of linear theory and a convolution on the line-of-sight with 
a velocity distribution. This model is accurate enough to robustly 
constrain the effect of massive neutrinos on RSD when applied on 
scales < 50 hT y Mpc . 

We have compared the findings from the ACDM and the 
ACDM+v simulations, a nd analysed their agreement with the ana - 
lytical predictions of ST dSheth et al.l200ll : lsheth & Tormenl2003) . 
Concerning the halo MF, we recover what theoretically expected, 



© 0000 RAS, MNRAS 000, 000-000 



Neutrino effects on LSS 9 




Figure 6. Two-point auto-correlation function, i~, in real and redshift-space. The contours represent lines of constant correlation, S,(r± 1 ni ) =0.15,0.25,0.5,1,3, 
for M v = (blue), M v = 0.3 eV (magenta) and M v =0.6 (red), respectively. Different panels show the results at redshifts z = 0.6, 1 , 2, as labeled. 



i.e. that, starting from the same A^(Ao) as initial condition, mas- 
sive neutrinos suppress the comoving number density of DM haloes 
by an amount that increases with the total neutrino mass M v . 
The suppression affects mainly haloes of mass 
10 M@/h, depending slightly on the redshift z. As an example, 
the number density of haloes with mass 1O 14 M /A at z = de- 
creases by ~ 15% for M v = 0.3 eV and by ~ 30% for M v = 0.6 
eV, and, at z = 1, by-~ 40% and ~ 70%, respectively. Moreover, 
with increasing z, the suppression of the halo number density due 
to free-streaming neutrinos moves towards masses M ^ 10 14 Mq//j. 

With regard to the halo clustering in the real-space, we observe 
that the trend of the halo correlation function t,(r) is opposite to the 
dark matter one. In fact, on one side, for a fixed A^»(&o) the total 
matter correlation function decreases with respect to the ACDM 
case due to neutrino free-streaming, in particular on small scales. 
On the other side, the halo correlation function undergoes the op- 
posite trend since the halo bias results to be significantly enhanced. 
For M v — 0.6 eV, we find that the halo correlation function in the 
presence of massive neutrinos at z = 1 is ~ 20% larger than in a 
pure ACDM model, and at z = 2 this difference rises up to ~ 40%. 
Also in this case, the theoretical ST bias model reproduces cor- 
rectly the numerical findings, inside the statistical errors, and, as in 
the ACDM cosmology, the halo bias results to be scale-independent 
on scales larger than r > 20 h Mpc . 

Considering RSD, we find that the rise of the spatial halo clus- 
tering due to massive neutrinos is less enhanced in the redshift- 



space than in the real-space. In fact, on large scales, the value as- 
sumed by the bulk flow, (v 2 (R t )), in a ACDM+v cosmology is 
smaller than in a pure ACDM one. On small scales, also FoG get 
decreased in the presence of massive neutrinos, so that the best-fit 
values of P and a 12 reduce by an amount which increases with M v 
and z. As an example, fixed the same initial condition on A^(fco), 
at z = 0.6 the (3 best-fit values decrease by ~ 10% for M v = 0.3 eV, 
and by ~ 25% for M v = 0.6 eV. Likewise, the O12 best-fit values 
decrease by ~ 25% for M v = 0.3 eV, and by ~ 45% for M v = 0.6. 

If not taken correctly into account, these effects could lead to 
a potentially fake signatures of modified gravity. Moreover, esti- 
mates of (3 and O12 can be used to extract measurements of the total 
neutrino mass and may help breaking degeneracies with the other 
cosmological parameters. 

However, these effects are nearly perfectly degenerate with the 
overall amplitude of the matter power spectrum as characterised by 
rj8. This strong M v -rJg degeneracy undermines the potentiality of 
the mentioned methods in constraining the neutrino mass. For in- 
stance, the difference between the halo MFs in the ACDM+v and 
ACDM models largely decreases if we normalise the two cosmolo- 
gies to the same rjg. Similarly, when analysing RSD, we find that 
the p best-fit values fall in the shaded grey bands of Fig. [7] rep- 
resenting the propagated ~ 10% theoretical bias error, and which 
contain the theoretical predictions obtained in a ACDM cosmology 
renormalised with the rjs value of the simulations with a massive 
neutrino component. For such a value of the bias error, we are pre- 
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Figure 7. Best-fit values of p-<Ji2, as a function of M v and redshift (points), compared with the ST theoretical predictions (solid lines). The dashed lines 
show the theoretical predictions in a ACDM cosmology with the same A%(ko). The blue dotted lines show instead the theoretical predictions for a ACDM 
cosmology normalised to the o~8 value of the simulation with a massive neutrino component, as explained in the text. The shaded grey bands represent the 
propagated ~ 10% theoretical bias error. The error bars represent th e scatter in the me asured (5 obtained dividing the simulation box in 27 sub-boxes, and 
rescaled by the square root of the total volume of the simulation box (Guzzo et al. 200^). 



vented to distinguish the effect of massive neutrinos on (3, if we use 
as initial condition the same rJs value both for the ACDM+v and 
ACDM cosmologies. 

Nonetheless, the rjg renormalisation of the matter power spec- 
trum does not totally cancel the neutrino effects which, in this case, 
depending on the M v value, alter the MF shape and amplitude es- 
pecially for the less massive objects. As an example, at z = the 
difference between the halo MFs with and without massive neutri- 
nos is ~ 3% at M = 10 13 M S jh for M v = 0.6 eV. The detection of 



this small effect depends on the sensitivity in measuring the halo 
MF at masses M < 1O 14 M jh. 

More promising are measurements of p. In Fig.[8]we show, as 
a function of M v and z, the relative difference between the theoret- 
ical P values calculated in the ACDM+v and ACDM cosmologies, 
normalised to the same rjg. At z = 1 and for M v > 0.6 eV, the rel- 
ative difference with respect to the M v = case is AP/P > 3%. 
This results is interesting, since future nearly all-sky spectroscopic 
galaxy surveys, like EUCLID, JEDI and WFIRST, should be able 
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Figure 8. The relative difference between the theoretical p values calculated 
in the ACDM+v and ACDM cosmologies, normalised to the same 08 . 

to measure the linear redshift-space distortion parameter with er- 
rors < 3% at z ^ 1, per redshift bin. This means that, even exploit- 
ing information from P measurements alone, they will contribute, 
along with other cosmological probes, to constrain the value of the 
total mass of cosmological neutrinos. 
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